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O^l ■ We study a coupled quintessence model in which the interaction with the dark matter sector is 

a function of the quintessence potential. Such a coupling can arise from a field dependent mass 
| term for the dark matter field. The dynamical analysis of a standard quintessence potential coupled 

Ifr . with the interaction explored here shows that the system possesses a late time accelerated attractor. 

In light of these results, we perform a fit to the most recent Supernovae la, Cosmic Microwave 
Background and Baryon Acoustic Oscillation data sets. Constraints arising from weak equivalence 
principle violation arguments are also discussed. 
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I. INTRODUCTION 



o 

i-C ■ Cosmological probes indicate that the universe we observe today possesses a flat geometry and a mass energy 
density made of ~ 30% baryonic plus cold dark matter and 70% dark energy, responsible for the late-time accelerated 



Q . expansion. Unveiling the origin and the nature of dark energy is one of the great challenges in theoretical cosmology. 
The simplest candidate for dark energy is the cosmological constant, which corresponds to a perfect fluid with an 
equation of state w = p/p = — 1. The ACDM model, i.e. a flat universe with a cosmological constant, is in very good 
agreement with current observational data. However, from the quantum field approach, the bare prediction for the 
\ current vacuum energy density is ~ 120 orders of magnitude larger than the measured value. This situation is the 
y—i • so-called cosmological constant problem. In addition, there is no proposal which explains naturally why the matter 
and the vacuum energy densities give similar contributions to the universe's energy budget at this moment in the 
cosmic history. This is the so-called why now problem. A possible way to alleviate this problem is to assume a time 
varying, dynamical fluid. The quintessence option consists on a cosmic scalar field cf> that changes with time and varies 
across space, and is slowly approaching its ground state. In principle, the quintessence field may couple to the other 
fields, see Refs. [H-ll4j]. In practice, observations strongly constrain the couplings to ordinary matter [l5|. However, 
interactions within the dark sector, i.e. between dark matter and dark energy, are still allowed. The presence of these 
interactions could significantly change the universe and the density perturbations evolution, the latter being the seeds 
for structure formation. We explore a scalar field dependent dark matter-dark energy coupling and confront the model 
predictions with current cosmological data. For models similar to the one studied here, see e.g. Refs. [Tl-fTTL Il3l. flij. 
The structure of the paper is as follows. Section |TT] presents the lagrangian theory responsible for the dark sector's 
coupling explored here. Section Hill describes the cosmological data sets used in the analysis. The dynamical stability 
of the coupled model and the fits to several cosmological observables are presented in Sec. HVl Weak Equivalence 
Principle violation constraints are explored in Sec. [Vj We draw our conclusions in Sec. IVI1 



II. COUPLED QUINTESSENCE 

Let us consider an interaction between the dark energy scalar field cf> and a cold dark matter field $ through the 
dark matter mass term md m (4>)^'^ ■ This form of interaction is inspired by the universal coupling to all species present 
in scalar-tensor theories in the Einstein frame [l[. Given that, observationally, interactions between the dark energy 
field and ordinary matter are strongly constrained [l5j , we will assume that the dark energy field 4> does not couple 
to baryons. At the level of the stress-energy tensor conservation equations, a dark energy-dark matter interaction 
m dm (<j))^>^ implies 

v,i? dm)t/ = pwfa^ = -v M ip e) „ with m = 9 J^M . (i) 
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The background evolution equations read 

$ + 3Hi + V,f = -P(<P)Pdm ; (2) 

Pdm+3Hp dm = P{4>)pdm4> , (3) 

where we have used a spatially-flat Friedmann Robertson Walker metric (ds 2 — —dt 2 + a 2 dx 2 ) and the dot refers to 
time derivative d/dt. 

Mostly all of t he p revious studies on coupled quintessence models have assumed that the coupling 0(<j)) is a constant, 
see e.g. Ref. 0, Il6| . In this paper, we consider a coupling varying with time (see also Ref. [T3, 13); given by some 
power of the potential of the dark energy field, i.e. /3(cj>) oc V(4>) n - For a slow-rolling quintessence field, the former 
assumption is equivalent to assume a coupling proportional to some power of the dark energy density /3(<p) oc p^ e . 
Therefore, our choice of interaction will naturally provide a dark sector interaction oc pdePdm if ti = 1. The model 
presented here should be understood as a lagrangian basis for more phcnomcnological approaches as, for instance, the 
one presented in Ref. [l4j . We illustrate the case of a (coupled) quintessence model, characterized by an exponential 
potential 

V(<j>) = M 4 exp{—a<f)/Mpi] with m dm = m exp (V (</>)/ p° cr ) n , (4) 

where p° cr is the current critical mass-energy density today, Hq = 87rG/3p° r = K 2 /3p° r , with G = \jM 2 v The scalar 
field dependence on the dark matter mass ensures that P(<j>) oc V(4>) n . 

For the stability analysis of our coupled dark matter-dark energy model with the potential of Eq. Q, we shall 
focus either on the matter-dominated era or on the late time dark energy domination period, neglecting the radiation 
contribution. Therefore, 

9 

Ky P 

ftdm + 0^ = 1 with fl = . (5) 
The baryons have been also neglected 1 . We introduce the dimensionless variables x, y, as in the uncoupled case (l9j | 



Jh 2 ' v ~ m 2 



x=T^> V=^H- (6) 



The positivity of the potential energy implies that y > 0. In the new variables x and y, the equations of state are 

2 2 

P<P x — V Ptot n 2 2 (~\ 

= — = — 5-, wtot = — — x - y . (7) 

P0 x l + y 1 p tot 

The condition for a late time accelerated expansion period is still u> tot < —1/3, as in the uncoupled case. The Hubble 
evolution equation can be written as 

|;=-5 (1 + * 2 -^- (8) 

The resulting evolution equations do not allow for a two-dimensional representation of this model since H can not be 
uniquely determined from the evolution Eqs. @ and ([3]) using exclusively the variables x and y. Following Ref [20| . 
we define a third dynamical variable z 

(9) 



H + Ho 

The condition < z < 1 ensures the compactness of the phase space. 

III. COSMOLOGICAL DATA USED IN THE ANALYSIS 

In this section we describe the cosmological data used in our numerical analysis. Three different geometrical probes 
(Supernovae la (SNIa), Cosmic Microwave Background (CMB) and Baryon Acoustic Oscillations (BAO) data sets) 
are exploited to derive the cosmological bounds on the coupled quintessence model. 



1 We neglect the presence of radiation and baryons in the stability analysis. For numerical purposes and fits to observational probes, we 
include both the radiation and the baryon contributions to the total mass-energy density. 
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A. The Supernova Union Compilation 

The Union Compilation 2 [21] consists of an update of the original Union compilation [22j with 557 SNIa after 
selection cuts. It includes the recent large samples of SNIa from the Supernova Legacy Survey and ESSENCE Survey, 
and the recently extended data set of distant supernovae observed with the Hubble Space Telescope (HST). In total 
the Union Compilation presents 557 values of distance moduli (/i) ranging from a redshift z of 0.015 up to z = 1.4. 
The distance moduli, i.e. the difference between apparent and absolute magnitude of the objects, is given by 

^ 51 ° g (^) +25 ' < 10 > 
where cIl{z) is the luminosity distance, di,(z) = c(l + z) J* H(z)~ 1 dz. The \ 2 function used in the analysis reads 

X 2 SNIai c i) ^^2(f^(Ct,z)~ ^obs(z))C;.l,(^(ci,z r ) - fJ, bs(z')) , (11) 

where, Cj refer to the free parameters of the coupled model and C is the covariance matrix with systematics included, 
see [2lf for details. 



B. CMB first acoustic peak 



We exploit the CMB shift parameter R, since it is the least model dependent quantity extracted from the CMB 
power spectrum [23| . i.e. it does not depend on the present value of the Hubble parameter H . The reduced distance 
R is written as 

/•1089 

R = (n m Ho) 1/2 / dz/H(z) . (12) 
Jo 

We use the CMB shift parameter value R — 1.7 ± 0.03, as derived in Ref. [23j, where it has been explicitly shown 
that the value of the shift parameter R is mostly independent of the assumptions made about dark energy. The x 2 
is defined as Xcmb( c i) = K R { c i) - Ro)/vr } 2 - 



C. BAOs 



Independent geometrical probes are BAO measurements. Acoustic oscillations in the photon-baryon plasma are 
imprinted in the matter distribution. These BAOs have been detected in the spatial distribution of galaxies by the 
SDSS J43 at a redshift z = 0.35 and the 2dF Galaxy Redshift Survey [H (2dFGRS) at a redshift z = 0.2. The 
oscillation pattern is characterized by a standard ruler, s, whose length is the distance that the sound can travel 
between the Big Bang and recombination and at which the correlation function of dark matter (and that of galaxies, 
clusters) should show a peak. While future BAO data is expected to provide independent measurements of the Hubble 
rate H(z) and of the angular diameter distance Da(z) = di,{z)/{l + z) at different redshifts, current BAO data does 
not allow to measure them separately, so they use the spherically correlated function 

/ cz \ 1/3 

Dv(z) = (^iOO^j) ■ ( 13 ) 

In the following, we shall focus on the SDSS BAO measurement. The SDSS team reports its BAO measurement in 
terms of the A parameter, 



A(z = 0.35) = D v (z = 0.35) ^""f" , (14) 

0.35c 



where A S dss(z = 0.35) = 0.469 ± 0.017. The X function is defined as 



x!ao( c *) = [{.Mpi,* = 0-35) - A SDSS (z = 0.35))/(j^ = o.35)] 2 - 
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IV. DYNAMICAL ANALYSIS AND COSMOLOGICAL CONSTRAINTS 

The coupling /3(4>) of Eq. ([T]), for the dark matter field dependence of Eq. ((4]), gives a dynamical coupling which 
reads 



(15) 



In the next sections, we will study the impact of such a coupling in the dynamical behavior of the system as well as 
in the numerical analysis performed to the cosmological data sets considered here. 

A. Stability analysis 

We study the dynamical behavior of the dark matter-dark energy dynamical system in the matter dominated period. 
In terms of the x,y and z variables defined in SecHH Eqs. ©, © and © read: 



x' = -3x + ~-£L y 2 + -4Ly 2 " (1 / )2 " (1 -x 2 - y 2 ) + -x(l + x 2 - y 2 ) 
4 4 z 2n V ' 2 V y ' 

y = ~ a 4^ xy + 2 y ^ +x ~ y ' ' 

z ='\z{l-z){l + x 2 -y 2 ). 



(16) 

(17) 
(18) 



The associated critical points are presented in Table HI They are independent of the value of n, assumed to be n > 1. 



Eigenvalues 



u>t Acceleration? Existence? 



3 3 _ 3 

2 ' 2 ' 2 



No 



V a 



±1 



3,3, 3t 



No 



Va 



4V / 37T 



^,^-3,sgn(-(-a 2 + 48^)")«) 



-1 + — 

- 1 ^ 247T 



Q < 167T 



< 48tt 



3 _ 3 3 
'2 ' 2 ' 2 



No 



V a 



±1 



-3, 3, 3=F ja 



No 



V a 



3 3(a+ v / -7ci 2 +192^) 3(-a + % /-7ct 2 + 1927r) 
2 ' 4q ' 4a 



24- 



No 



a 2 > 24n 



4V^ 



16 - !f 

4 V 37T 



" 16tt ' ° "•" 16ir ' ° "T" 8?r 



-1 + — 

- 1 ^ 247T 



Q < 167T 



< 48tt 



TABLE I: Critical points and associated eigenvalues for the exponential potential model of Eq. Q and n > 1. 



Notice that our results are very similar to the ones obtained in model C of Ref. 14| 2 . The main difference among the 
results presented here and those presented in Ref. [14[ lies in the eigenvalues for the z = case. The critical points 
for a matter dominated period followed by an accelerated expansion are however rather equivalent. 

The exponential coupled model of Eq. (U) allows for a matter dominated era at early times, corresponding to our 
first critical point, see Tab. HI This critical point is an unstable fixed point regardless of the value of the exponential 
potential parameter a. The last critical point of Tab. Q]is an accelerated attractor for a 2 < 16tt. We show the phase 
space trajectories pointing towards this attractor in Fig. [1] for z — 1 and a = 1. 



2 In Ref. [14|| . the coupling was chosen to be proportional to the Hubble rate parameter Hq. 
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FIG. 1: Phase space trajectories for the exponential potential studied in this paper at z — 1. The plot illustrates the stable 
node located at (x,y) — (0.16,0.99) for a = 1. This stable node corresponds to an accelerated attractor and its existence is 
independent of the value of the power n of the potential which appears in the dark sector interaction. 



B. Cosmological constraints 

We have already shown that the background dynamics of the exponential potential coupled model studied here offers 
a suitable framework to describe the late time accelerated expansion of the Universe. We now present the constraints 
which arise from the data sets described in Sec. IIIII Both baryon and radiation contributions to the expansion rate 
have been included in the following analysis. In the discussion, we make use of the individual chi-square functions, 
and the global chi-square is defined by 

Xtot(Ci) = XSNIa( c i) + X 2 BAo( c i) + XCMs( c i) ) (19) 

where a refers to the free parameters of the coupled model under study. The coupled model analyzed here contains 
three parameters a,n and M. Two parameters determine the scalar field exponential potential: its amplitude is set 
by the mass scale M and a appears in the argument of the exponential. The parameter n fixes the power of the 
scalar field potential appearing in the dark sector interaction. Notice that, in this interacting quintessence model, a 
multiplied by n plays the role of a dimensionless coupling which sets the magnitude of the interaction in the evolution 
equations, see Eq. (fT5j). 

Figure [2 left panel, shows the results of the x 2 analysis to SNIa data. In the right panel, the results for the global 
fit analysis are depicted. From left to right, the curves show the 68.3 and 95.4% C.L. allowed regions for the n = 1, 
n = 2, n = 4 and the uncoupled (f3((j>) = 0) cases. For this particular analysis the initial conditions for the scalar 
field are set to cf>i n — M p i and <f>i n — 0. We have checked the robustness of our results versus the scalar field initial 
conditions. Our conclusions remained unchanged, in agreement with the results of Ref. [26| . 

In general, larger values of a imply a smaller scalar potential. To compensate this effect, larger values of the ampli- 
tude of the potential M are needed. This explains the shape of the degeneracy between a and M in Fig. [2], being these 
two parameters positively correlated. For small values of both M and a, for instance, (M, a) ~ (1.8 10 _31 M p /, 0.1), 
the ratio V/p® r appearing in the coupling term is smaller than one. When the value of n is increased from 1 to 4, 
the uncoupled case behavior is recovered due the suppression of the coupling term (V/ p a cr ) n . Indeed, notice that the 
n = 4 curve is almost superimposed to the uncoupled quintessence curve in Fig. H in the low (M, a) region. 

For larger values of a and M (within the allowed regions), the ratio V/ p° cr increases. The dark sector interaction 
becomes the dominant source term for the scalar field evolution in Eq. @. In addition, the dark matter energy 
density, proportional to exp [(V \(j>) j p° r ) ] (see Eq. Q), starts to dominate the total energy density. Such a large 
contribution from the dark matter energy density is not compatible with SNIa data. This is precisely the reason for 
the bending of the curves in the left panel of Fig. [2] Notice that the shape of the curves for the coupled cases differs 
significantly from the uncoupled case for larger values of a and M (within the allowed regions). The turn- over of 
the coupled model curves occurs at different values of the parameters M and a, depending on the value of the n 
parameter. Namely, for n = 4 the turn over shows up at smaller values of M and a than those corresponding to the 
n = 1 case. This is due to the fact that the /3(</>) term given by Eq. @ is proportional to n and therefore the strength 
of the coupling term in this region of the M and a parameters grows with n. 

Notice, from the global fit results of Fig. [5] (right panel), that the allowed regions for large values of n, a and M 
become significantly smaller than those arising from a fit exclusively to SNIa data. This is due to the CMB constraint 
which tends to favor smaller values of M in both uncoupled and coupled cases. Indeed, larger values of M are 
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FIG. 2: (Left panel) Analysis of the coupled exponential potential model of Eq. Q). The contours denote the 68.3 and 95.4 % 
C.L. allowed regions arising from a fit to SNIa data. The first three curves from left to right depict the results for n — 1,2 and 
n — 4. The last curve depicts the results arising from the fit to the uncoupled model. (Right panel) Same as in the left panel, 
but using the global fit results, i.e. the combined analysis to SNIa, CMB and BAO data sets. 



associated to very small values of the dark matter energy density, which directly influences the CMB shift parameter, 
see Eq. fT2|), 

The best fit point for the uncoupled case is located at (M,a) — (2 10 _31 M p ;, 0.25) and it corresponds to Vide = 
0.72 and Qdm = 0.28, being fi^e and Qdm the current values of the dark matter and dark energy energy-densities, 
respectively. The equation of state of the dark energy field is = —1. The associated global \ 2 is 562 compared 
to the 557 effective degrees of freedom in our analysis. While for the n = 1 coupled case the fit becomes weaker, it 
improves for larger values of n. For the n — 4 coupled scenario, the best fit is associated to a x = 542. However, 
the best fit point for the n = 4 case corresponds to a cosmological constant scenario, since it is located at a — and 
M = 1.85 10~ 31 M P 2, with the derived cosmological parameter values fide = 0.67, ftdm — 0.33 and = —1. Indeed, 
for ACDM universe we obtain x 2 = 531. 

In summary, the cosmological data sets considered in the current analysis favor the large n regime in a region of 
the (M, a) plane where both the coupling term and the dynamics of scalar field potential are negligible. In the next 
section we will explore the constraints from weak equivalence principle violation arguments. 



V. WEAK EQUIVALENCE PRINCIPLE CONSTRAINTS 

Let us consider an interaction between fermionic dark matter, \f r , and a light pseudo scalar boson, cj>, that interacts 
with the dark matter through a Yukawa coupling with strength g, described by the lagrangian 



L = i^V^* - m^* - -V^VV - V{4>) + g<j>^^ , 



(20) 



where is the dark matter mass (independent of the scalar field). For g ^ 0, on scales smaller than r s — m^ 1 , the 
Yukawa interaction acts like a long-range 'fifth' force in addition to gravity. The effective potential felt between two 
dark matter particles is 



V(r) = - 



1 + aYuk exp 



(21) 



with 



"Yuk 



4tt m 2 



(22) 
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The authors of Ref. [27| studied the impact of such a long-range interaction on both galaxy cluster masses and dark 
matter growth, reporting an upper bound of ayuk < 1-3 (for oy u k > 0). 

Kesden and Kamionkowski (K&K in the following) [H, [HI , analyzed the consequences of Weak Equivalence Princi- 
ple (WEP) violation for dark-matter on galactic scales, focusing on dark-matter dominated satellite galaxies orbiting 
much larger host galaxies. They concluded that models in which the difference among dark matter and baryonic accel- 
erations is larger than 10% are severely disfavoured. The relevant parameter constrained in K&K analysis is y / OYuk- 
For reasonable models of the Sagittarius satellite galaxy tidal stream, K&K found an upper bound of oy u k < 0.04, 
corresponding to 



g/m^ < 4.2 x 10~ 20 GeV" 



(23) 



Since the K&K limit turns out to be the most stringent one, we exploit it here to set constraints on the exponential 
coupled model explored along the paper. We expand the interaction term m c ^ m (</>) , ^ , ^ , , keeping only the linear terms 
in the scalar field. The former approximation is possible due to the fact that acj)/M p i <C 1 within the viable regions 
determined in the previous section. Therefore the linear term will clearly dominate the dark matter-scalar field 
interaction. We can write Eq. ^ as: 



U s v> 1 (V{<j>)\ 



k=0 



V{4>) 



M 1 - a-. 



M, 



pi 



(24) 
(25) 



It is straightforward to obtain an expression for both the fermion mass m$ and the coupling g in the Yukawa interaction 
term as defined in Eq. ([2"0"]l 



mo exp 



M 4 

Per 



na ( M 



M p i \ p cr 

Therefore, if we express the amplitude of the potential as M = Ax 10~ 31 M p i, the limit of Eq. (f23| becomes 



a < 



0.1a 



>&TT A 4 > 

v 3 100/ 



(26) 
(27) 

(28) 



For A = 1.85, which lies in the range of the viable regions obtained in Fig. [21 the former bound translates into a 
bound on the parameter a of the exponential coupled model of a < 0.51, 0.26 and 0.14 for n — 1, 2 and 4, respectively. 
These bounds are stronger than those arising from cosmological observations and further restrict the allowed regions 
shown in Fig. [5J In addition, WEP bounds are complementary to cosmological constraints, since they have opposite 
trends: while cosmological bounds are rather loose when the amplitude of the potential increases, WEP limits get 
much stronger. Notice however that WEP bounds have been obtained using the strongest fifth force constraint, i.e. 
using the K&K limit. Mildest bounds on coupled quintessence models will arise if more conservative WEP bounds 
are applied. 



VI. CONCLUSIONS 



We have studied a time varying-interaction among the dark matter and the dark energy sectors. The non-minimally 
coupled dark energy component is identified to a dynamical quintessence field, and it is coupled to the dark matter 
field via the dark matter mass term. The form of the interaction has been chosen to ensure an energy exchange 
between dark matter and dark energy proportional to the product of the dark matter energy density and the nth 
power of the scalar field potential. For a slowly rolling scalar field and n = 1, the model presented here provides 
a possible effective lagrangian description of pure phenomenological quadratic interacting models such as the one 
studied in Ref. [HJ. 

The form for the scalar field self-interacting potential is assumed to be an exponential function of this field. The 
model has then been shown to possess a late time stable accelerated attractor regardless of the value of the n parameter. 
We have also explored the constraints on this interacting model arising from the most recent SNIa, CMB and BAO 
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data. While the fit improves slightly when allowing for a dark matter-dark energy interaction, the best fit point lies 
in a region in which both the coupling and the dynamics of the field arc negligible. Therefore, current data does not 
favor this coupled dynamical model, even if it involves more parameters than the simplest cosmological scenario, i.e. 
a ACDM universe. A coupling between the two dark sectors can also be constrained by Weak Equivalence Principle 
violation arguments. We have derived the constraints arising from fifth-force searches. For the exponential potential 
coupled model studied in this paper, WEP constraints are very strong and much tighter than cosmological bounds. 
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